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Abstract: The auxiliary /dynamic decoupling method of hep-th/0609001 applies to per- 
turbations of any co-homogeneity 1 background (such as a spherically symmetric space-time 
or a homogeneous cosmology). Here it is applied to compute the perturbations around a 
Schwarzschild black hole in an arbitrary dimension. The method provides a clear insight 
for the existence of master equations. The computation is straightforward, coincides with 
previous results of Regge- Wheeler, Zerilli and Kodama-Ishibashi but does not require any 
ingenuity in either the definition of variables or in fixing the gauge. We note that the 
method's emergent master fields are canonically conjugate to the standard ones. In addi- 
tion, our action approach yields the auxiliary sectors. 
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1. Introduction and Summary 

Recently an improved understanding of perturbations around backgrounds with (at most) 
one non-homogeneous coordinate ( "co- homogeneity 1") was achieved 0. Two cases of 
physical interest fall under the title of "co-homogeneity 1 backgrounds" : static spherically- 
symmetric space-times and homogeneous cosmologies. In this paper we apply the method to 
the Schwarzschild black hole in an arbitrary dimension, re-deriving the master equations for 
the perturbation spectrum and wave profiles while obtaining new insights to be mentioned 
below. 

The method of Q includes the method known as "gauge invariant perturbation theory" 
together with an explicit treatment of the auxiliary sector hence we refer to it as the 
" auxiliary / 'dynamic decoupling method". In addition it was proven in greater generality, 
demonstrating that the essential condition is the dimension of co-homogeneity (Id). 

Before discussing the specifics of black hole perturbations let us recall the basic ingredi- 
ents of the general theorem |TJ . By performing dimensional reduction over the homogeneous 
coordinates we may consider the problem to be essentially Id. It is found that a Id ac- 
tion with uf fields invariant under no gauge functions can be decoupled into two sectors, 
auxiliary and dynamic, where the auxiliary sector is purely algebraic as its name suggests 
and all fields in both sectors are gauge invariant. The auxiliary sector consists of ng fields 
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which are nothing but the (shifted) metric and vector fields from the Id perspective, and 
hence both have formally minus one degrees of freedom. The dynamic sector contains 
hf — 2 uq fields which remain after the gauge functions are eliminated from the np — uq 
non-auxiliary fields. All field redefinitions are local in Id and invertible. 

In addition to the decoupling between auxiliary and dynamic fields there is another 
useful decoupling mechanism. From the Id perspective the isometry group becomes the 
theory's gauge group under which the fields are charged. Therefore, the quadratic action 
decouples according to representations of this group. 

Chronologically, the general theorem of followed by generalizing "the" derivation 
of the Gregory-Laflamme mode [§] (see [||, U for reviews of the black-hole black-string 
transition as well as ||). 



Let us return to discuss the specific case of perturbations around Schwarzschild. Of 
course, the results for this case were already known. The 4d case is by now classic: the odd 
(vector) sector reduces to the Regge- Wheeler master equation || while the even (scalar) 
sector reduces to the somewhat more cumbersome Zerilli master equation 0]. The higher 
dimensional generalization however, is relatively recent, and was achieved by Kodama 
and Ishibashi [|| (see also || for the master equation of static perturbations through a 
different method, and |[(]] for an application to a stability analysis of higher dimensional 
Schwarzschild ). 

Our method reproduces and confirms these results while adding the form of the auxil- 
iary sectors and providing the advantages of a different point of view (somewhat in analogy 
to Moncrief's method in 4d, see also [12, ^] for some recent discussions of perturba- 



tions of the Schwarzschild black hole in the gauge-invariant formalism). The method has 
the advantage of being practically straightforward. We use an action approach and sim- 
plify it in a natural way. We explain why the equations in each sector can be reduced to 
a single master equation (as anticipated in |l]]). Moreover, there is no need for ingenious 
field redefinitions as in [|| nor is there a need for choosing a gauge-fix as in || . In addition 
our treatment is independent of the choice of coordinates for the background Schwarzschild 
space-time. 

While our method identifies the space of dynamic fields, it only fixes the definition 
of the dynamic fields up to a linear 1 , radially-dependent canonical transformation. It is 
found that the action is considerably simpler for some choices, and the standard master 
fields of Regge- Wheeler, Zerilli and Kodama-Ishibashi are essentially canonically conjugate 
to the ones which naturally emerge from our method. The standard master fields have the 
additional property that their r and t derivatives combine in the master equation to the 
form of a Laplacian in the reduced 2d (r, t) space. This property comes as a surprise from 
the current Id perspective which concentrates on the r coordinate, and to gain insight into 
it, it is probably necessary to study the properties of perturbations on co-homogeneity 2 
spaces. Another open question from our perspective is whether the standard form of the 
equations is in some sense optimal or whether further simplification is possible, especially 



To preserve the linearity of the equations of motion. 
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in the scalar potential which is quite involved. 

We start in section [2] with a review of the auxiliary/dynamic decoupling method. In 
section ^ we derive the quadratic Lagrangian for the perturbations. While the derivation 
of this Lagrangian is completely straightforward, it is quite lengthy, and perhaps there is 
room to make it more economical. Due to the decoupling of representations at quadratic 
order (mentioned above) the Lagrangian decouples to three sectors known as tensor, vector 
and scalar, and which reduce in the 4d case to the even and odd sectors (the tensor sector 
degenerates) . 

In section |] we proceed to apply the method to this Lagrangian. The tensor sector 
has one field and no gauge functions (rip = 1 uq = 0) and hence requires no treatment. 
The vector sector has np = 3 uq = 1 and since 3 — 2*1 = lwe immediately expect a 
single dynamic field, and hence a decoupled equation of motion, namely a master equation. 
Indeed our straightforward manipulation results in the dynamic and auxiliary sectors of 
the action. The dynamic sector encodes a master equation which is shown to be equivalent 
to the Kodama-Ishibashi equation (or Regge- Wheeler in 4d). The scalar sector is more 
involved: it has np = 7 no = 3 but as 7 — 2*3 = 1 there is still a single master equation. 
Again we follow our method to derive the auxiliary sector together with a master equation 
and the latter is confirmed to be equivalent to the Kodama-Ishibashi equation (Zerilli's in 
4d). Our main results are summarized in section [f| 

It would be interesting to apply the auxiliary/dynamic decoupling method to other 



space-times including perturbations around the critical merger solution [14]. 
2. Review of the auxiliary /dynamic decoupling method 

Consider perturbations around a space-time X which is of co-homogeneity 1, namely that 
all but one of its coordinates (denoted by r) 2 belong to a homogeneous space. One starts by 
effectively reducing the problem to a Id problem on r as follows. The metric perturbation 
tensor hav decomposes into several fields according to whether fj, is r or belongs to the 
homogeneous space. Furthermore, each component of h^ u can be expanded into harmonic 
functions of the homogeneous space in order to achieve a separation of the variables. For 
instance, for the time translation symmetry U(l)t we expand in exp(itvt), while for the 
spherical symmetry SO(D — 1)q we expand in spherical harmonics (of the appropriate 
tensor type). 

Altogether the action reduces to a Id action in r, with a certain field content 

= <^(r) i=l,...,n F (2.1) 

where np is the number of (real) fields. Actually the field content necessarily includes the 
Id metric h rr as well as massless vectors which originate from the isometry group of X and 
become the gauge symmetry in Id. Hence the dimensionally reduced theory is actually aid 



2 In |l] the essential coordinate was denoted by x, but here since we specialize to the spherically symmetric 
case we denote it by r. Some additional notation which differs between this paper and Jlj is: "DG" which 
stood for "derivatively-coupled" sector was exchanged for "A" which stands for "auxiliary" or "algebraic" , 
and "X" which stood for the dynamic sector was traded with "Dyn" . 
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gauged- gravity theory with additional matter content (which generically includes r-scalars 
as well as massive r- vectors and r-tensors). 

The invariance of X under infinitesimal diffeomorphisms £^ translates into invariance 
of the Id action under certain gauge functions, whose number we denote by n G , as follows 

t = e(r) a = l,...,n G (2.2) 

where (Gi(r), Go(r)) are a pair np x n G r-dependent matrices. 

The action is quadratic in <j> since we are interested only in linear perturbations, and 
we furthermore assume it to be quadratic in derivatives as is the case for classical GR. 

The standard counting of equations proceeds as follows. One uses the n G gauge func- 
tions to set n G gauge-fixings relations over the fields, possibly eliminating n G fields. The 
gauge-fixing is accompanied by n G constraint equations which are usually of lower order in 
derivatives as well as rif — n G other equations. |jj proved that under generic conditions 3 it 
is always possible in Id to completely eliminate the gauge through a transformation of the 
fields which is local (in r) and invertible, and that the action decouples into an auxiliary 
sector consisting of n G fields with an algebraic action and a dynamic sector consisting of 

n F -2n G (2.3) 

fields (generically second order in derivatives). Comparing to the standard counting we 
find that in Id it is always possible to have algebraic equations for both the n G gauge-fixing 
relations as well for n G of the equations of motion. 

As the full proof is given in [|J , we choose to limit ourselves to a description of the key 
mechanisms at work. The essential idea is a sort of "derivative counting" applied to the 
equation which states that the equations of motion (EOM) are satisfied for a pure-gauge 
field configuration. We define "derivatively gauged" fields or "A-fields" 2 to be the image of 



G\ (2^2). Physically, these are nothing but the r-tensors and r-vectors. The term with the 
highest number of derivatives (3) originates from the second-derivative terms in the EOM 
applied to the A-fields. Since this is the only term at that order we find that it must vanish, 
namely that the A-fields do not appear in action terms with two derivatives. Moreover, by 
considering a similar expansion of the second variation of the action we find that action 
terms in which only A-field appear are purely algebraic thus justifying the name "auxiliary 
sector". The relevant algebraic matrix is denoted by La- In the generic case that La is 
non-degenerate we may "complete the square" thereby decoupling the A-fields from the 
rest of the action after a suitable shift in their definition. 

Then we inspect the action of the gauge transformation on the new fields. The algebraic 
nature of the action in the auxiliary sector means that no symmetry is possible in that 
sector and hence no gauge invariance. By definition, the non-auxiliary fields are gauged 
only algebraically (£' does not appear in 6(j>) which means that the gauge can be completely 



3 Namely, that Gi is non-degenerate, as well as Gx and Lug = La which are respectively the gauge 
transformation when restricted to the non-auxiliary sector and the (algebraic) action of the auxiliary sector, 



both as defined in GJ. 
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eliminated by an algebraic re-definition of fields. Thereby np — Iuq fields are present in 
the other, dynamic sector, and completing the general description of the method. 



3. The Id quadratic action 

We consider gravitational waves in the Schwarzschild background in an arbitrary space- 
time dimension D. The isometry group of this space-time is 0(D — 1)q x U(l)t X ^2,T, 
where 0(D — 1)q is the spherical symmetry, U(l)t is time translation and ^2,T is time 
reversal. These isometries allow the use of the following "maximally general ansatz", 
namely an ansatz such that all of Einstein's equations can be derived from the variation of 
the gravitational action with respect to its fields 

ds 2 = e 2a ^ dt 2 + e 2b ^ dr 2 + e 2c « <m 2 D _ 2 , (3.1) 

where dQ 2 - ) _ 2 is the standard metric on the D — 2 sphere. The standard Schwarzschild 
coordinates are 

e 2a = f(r) = 1 - (r /r) D - 3 , e 2b = f~\ e 2c = r 2 , (3.2) 

but we may consider other coordinate systems (namely the choice of r) such as in ||15U , and 
for generality we shall not restrict to any specific choice. A general transformation from 
Schwarzschild coordinates shows that for any choice of coordinate we have 

c , = e a +b ^ c ^ a , = _ (jD _ 3)e 6-c sinha ; (3 3) 

where prime means a derivative w.r.t. r. 

The metric in the gravitational wave is = gffi + h^, where gffj stands for the 
background metric and h^ v is a perturbation. In order to write down the action that 
describes gravitational waves in a given background to linear order it is convenient to 
define the following tensors: 

V := V - \a^h\ , h:=Hfc, (3.4) 

and 

7V : = + V » h l ~ V " V) » (3-5) 

which is similar to the linearized part of the Christoffel symbol only with h — > h. Using 
these definitions the quadratic action can be compactly written as 

S = J d D x^ (2 j\ v >T X + -±—V^H) . (3.6) 

The gauge invariance of the full Einstein-Hilbert action has two consequences for the 
quadratic action: it is invariant both under a change of background coordinates (since it 
is manifestly invariant under local Lorentz transformations) and under the following gauge 
symmetry acting on the fields h^ v 

<*V = + V ^m> ( 3 - 7 ) 
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where is an arbitrary vector field. These two symmetries can be understood in a fol- 
lowing way: as the components of the metric in the Einstein-Hilbert action are expanded 
around the background, also the gauge functions should be expanded; an invariance under 
reparameterizations appears in the zeroth order and the invariance (3.7) appears in the 
first order. 

In this section we outline the derivation of an action in the radial direction, via a 
dimensional reduction of ( |3.6| ), and the details of the derivation are presented in Appendix 
|A|. We distinguish the components of the metric perturbation h^ v according to their tensor 
type, both with respect to r and with respect to the angular coordinates. Namely, in the 
radial direction h rr is an r-tensor, h r t and h r i (we denote by Latin indices the angular 
coordinates) are r-vectors, and the rest are r-scalars. Similarly on the sphere htt, h% r and 
h rr are fi-scalar-fields, ht% and h r i are fi-vectors-fields and hij form an O-tensor-field. 

As usual during dimensional reduction we expand each field into harmonic functions 
with respect to the isometry group. Actually the isometry group becomes the gauge group 
from the Id (radial) point of view, and each field is labeled by its representation under 
this group. For the spherical coordinates we need the scalar, vector and symmetric tensor 
spherical harmonics. The properties of these harmonics and the conventions used here are 
discussed in appendix ^. As discussed there, a spherical harmonic of a given type is labeled 
by a non-negative integer number I that determines the Dynkin indices [Zi,...Z n ] OI * the 
irreducible representation of SO(D — 1), and individual functions within the representation 
are further labeled by mi, ...m n which we denote collectively by m. There is a single family 
of scalar spherical harmonics, namely belonging to the representation [Z, 0...0], two 
families of vector spherical harmonics, namely Yf lm in the representation [7,0... 0] and 
YY lm in [Z, 1, 0...0], and 4 families of tensor spherical harmonics: Y$ lm and Y^ lm in 
[Z,'o...O], YY lm in [1, 1.0...0] and Y? lm in [Z,2,0...0]. The superscripts 5, V, and T stand 
respectively for scalar, vector and tensor and denote the three relevant families of irreps 
[Z, 0...0], [Z, 1, 0...0] and [Z, 2, 0...0], according to the lowest spin field for which the specific 
irrep family appears. Thus we should distinguish between two angular tensor types. The 
first tensor-type is with respect to local diffeomorpshims under GL(D — 1) is denoted, 
as was mentioned above, by J7-scalar -field, Q-vector-field etc. The second tensor type is 
with respect to the SO(D — 1) global isometries and will be denoted by S, V and T as 
above. From hereon we shall mostly refer to the latter tensor type (which specifies the 
representations of fields under the gauge group from the Id perspective), and when there 
is a possibility of confusion we shall call such types SO(D — l)-scalars, vectors or tensors. 

Dimensional reduction with respect to t means a Fourier expansion. Altogether the 
expansion is 

htt = diO ' l fy,mi...m„^i...m„ e! ' Jt ) 



l,m\...m n 



h^tf — ^tv ^ ' hjyp — hyty, ~Y c 



{hf Yf + hf Yy) e«" , h ri = {hf Yf + hX YY) e 
^ + h s Y^ + h v YY + h T Yf J 



h i3 = (h s Yfj + h s Yfj + h v YY + h T Y?) e iujt . (3.8) 
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where only the first equation shows explicitly the summation over representations I, over 
the indices m = m\...m n and the oj integration, and it should be understood to be implicit 
in all the other equations as well as those to follow. 

The notation for the coefficients here is the following: The superscripts denote the 
SO(D — 1) tensor-type (or irrep family) of the corresponding spherical harmonic, and the 
subscripts give the non-angular indices of h^. 

The gauge variations of h^ iu (3/7) allow one to calculate the gauge variations of the 
coefficients in these expansions if one expands also the components of the vector £^ into 
the spherical harmonics and Fourier expands in time: 

Zt = S?Y s e iut , tr = t S r Y s e^\ & = (£ s Y? + £ v Y?) e iut , (3.9) 



where the notation is similar to Q3.8| ). We shall spell out the gauge variation of the fields 
after dimensional reduction in the next section. 

The next step is to plug the expansions ([T^) into the action and calculate the angular 
integrals (all necessary formulae are given in Appendix ^|). In the next section we present 
the resulting action sector by sector in ( |4.1| , 4.4 , fO| , 4.31| ) and proceed to further analyze 
and decouple it. 



4d 

Our above-mentioned considerations do not apply directly to 4d since the rank of 0(3), the 
rotational symmetry group, is smaller than 2 and the above field-theoretic definitions of the 
vector and tensor families become meaningless since they are defined through the Dynkin 
indices [1, 1, . . . 0] and [1, 2, . . . 0]. However, since the decoupling into sectors relies mostly 
on the orthogonality relations of the spherical harmonics, the decoupling continues to hold 
in 4d, though with some changes: the scalar sector becomes the even sector, the vector 
sector becomes the odd sector and the tensor sector disappears, as we now proceed explain. 
Moreover, we find that the final form of the decoupled action in 4d can be obtained simply 
by setting D = 4 in the expressions we have for an arbitrary D. 

One change is that in 4d there is no tensor sector. To understand that we note that 
on S 2 the metric has only 3 components, and therefore the general decomposition of the 
metric into 4 components has to degenerate, and it turns out that the tensor component 
is the one that vanishes identically. 

The scalar and vector sectors still exist, but their field-theoretic definition is different. 
The group of isometries 0(3) includes also a ^2 of reflections in its center. Therefore 
any representation can be classified as either even or odd, and the quadratic action must 
decouple into these two sectors. Actually we assign a parity to each field, according to the 
parity of its I = component, and for higher I the parity gets multiplied by (— ) l . In this 
way the fields hf t , hf r , hf r , h^ , fyp , h s and h s are found to have even parity, whereas h\ , 
and h v are odd. Equivalently, the scalar sector becomes the even sector in 4d and the 
vector sector becomes the odd sector. 

The case of 5d is also special from the group theoretical point of view since SO{4) de- 
composes into SU(2) x SU(2). However, in this case no sector degenerates, and our general 
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theory continues to hold, since the rank is 2, and due to the reliance on the orthogonality 
relations. 



4. Decoupling each sector 

The quadratic action decouples into three parts: scalar, vector and tensor 

S = S S + S V + S T , (4.1) 

due to the orthogonality property of the spherical harmonics. In this section we shall 
further decouple each one of the sectors separately. 

In each sector the general formula for counting dynamic fields (2.3) leads us to expect 



a single dynamic field, which is called "the master field". This master field is unique up 
to an r-dependent factor, which we attempt to choose so as to simplify the action, but 
we do not pretend to optimize this choice here. Moreover, one is allowed to use canonical 
transformations which are non-Lagrangian (namely, the new field is a function not only 
of the old field, but also of its momentum). And indeed, we find that in order to obtain 
the standard form of the master equations we need to transform to a new field which is 
essentially the momentum conjugate to our original field. 

Another general feature is the r-tensor-type. All master fields are r-scalars, while 
the auxiliary fields are r-vectors and r-tensors. As usual, local Lorentz invariance in Id 
determines the required power of g rr = e~ 2b for each term, so we should be able to present 
the action such that b does not show up explicitly. To this end we use the following standard 
notation for derivatives and the integration measure 

d r = g rr d r = e~ 2b d r ds := e b dr . (4.2) 

Finally, all fields are of a definite length dimension. The dimensions of r,co are 1,-1 
respectively, while the Lagrangian has dimension —2 as a curvature. The dimensions of 
the metric functions and fields are 
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One can easily check that our Lagrangian summarized in (|4. 
dimension —2. 



and (|4.3l[) is of 



4.1 Tensor sector 

There is only one field in the tensor sector, namely h T . Naturally, it is gauge-invariant. 
The action simplifies by the following field redefinition 



~ 2c h T 



(4.3) 
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which is nothing but the divergenceless part of the mixed index components . In terms 
of <f> T the tensor action is given by 

S T = J e a+ ^ D -^ c dsC T 

C T = -I(5 r ^a> T + (C 2 e- 2c + W 2 e- 2a )0 T2 ) , (4.4) 

where 

C 2 :=l(l + D - 3) (4.5) 

is a value of a quadratic Casimir operator in a scalar representation with highest weight I 
(for more details about various representations see Appendix [(]). We will also find it useful 
to define 

C 2 :=C 2 -(D-2) = (l + D-2)(l-l). (4.6) 



In the expression (4-.L ) and throughout the rest of the paper bilinear expressions such as (f) 2 
must be understood as <f><j)* and the real part extracted. (Even if the two fields are different, 
the expression is the same for either possibility of which field to conjugate). 

4.2 Vector sector 

In the vector sector there are 3 fields: \i( , h% an d h . 4 The vector part of the action is 

S v = C 2 Je a+ ^ c dsC v 

C v = C 2 + d + C , (4.7) 

where C 2 is a "kinetic term", C\ is a term with one derivative w.r.t. r, and Cq is a 
"potential" : 

C 2 = - e ~ 2a ~ 2c d r hY 8 r hY - ^e" 4c d r h v 8 r h v , (4.8) 
A = -2iuo e - 2a ~ 2c d r h v t h v r + C 2 e' Ac d r h v h Y r , (4.9) 

£ = -e~ 2c (C 2 e~ 2c + ui 2 e~ 2a ) h v r h V r + Aiu> e'^ ^ h v r h\ 

-2C 2 e a - hc ^Th v r h v - e~ 2a - 4c (C 2 + 2{D - 3)e 2a ) hf 

2 



+ iC 2 u;e- 2a - 4c h v t h v - 9l e -*asc ^2 _ 4g 4a + 2{D _ 3 ya^ h v* (410) 

The gauge variations of the fields are (see ( |3,9| ) for the definitions) 

ShY = iuZ v , 5h? = Z v ' - 2c'i v , 6h v = 2i v . (4.11) 



4 In this subsection all fields will have a V superscript. It would have been possible to omit it, but we 
didn't do so in order to avoid a possible confusion with scalar fields in the next subsection with an otherwise 
identical notation. The same comment applies to the S superscript in the next subsection. 
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Altogether the number of fields and gauge functions is rip = 3, no = 1 and hence according 
to the general formula ( |2.3| ) we expect a single dynamic field (3 — 2*1 = 1) and one auxiliary 
field. 

since h v r is an r-vector it is seen that its variation includes £ and hence it is 
"derivatively-gauged" in the terminology of M. According to the general theorem for- 
mulated and proved there and reviewed in section [2| the derivatively gauged fields decouple 
from the rest of the fields after a certain shift in their definitions (completion of a square) 
and the newly defined fields are auxiliary. An important role is played by the quantity in 
the following definition 

A v := C 2 {1) e~ 2c M + to 2 e~ 2a ^ , (4.12) 



which comes from the auxiliary action, and in the scalar sector will be seen to generalize 
to be essentially the determinant of the auxiliary action (up to a multiplicative factor). It 
appeared in Q4.1C| ) as a coefficient of 2 ■ Now the shifted field can be written as 



H 



v 



-A v h^ + iive- a (2e 



e~ a dr) h 



v 



:Coe 



a-2c 



(2e- 



e- a d r ) h v 
(4.13) 



and derivatives of this new field do not appear in the Lagrangian, so this field is indeed 
auxiliary. 

The remaining part of the Lagrangian contains only algebraically gauged fields (fields 
whose variation does not contain £ Vr ). Hence, in accordance with the general procedure 
of each gauge function can be used to eliminate one field ("each gauge function shoots 
twice" - the gauge function already "produced" an auxiliary field). 

In our case we have single gauge function £ , so we can eliminate one algebraically 
gauged field out of two, or in other words, we can write the Lagrangian in terms of a single 
gauge invariant combination of fields. We define this dynamical field to be 



lV 



b\ = e~ c (u h v + 2ihY) 



(4.14) 



The Lagrangian then includes terms proportional to (j)' 2 , <j><fi', <j) 2 and 2 . The term <jxp' 
can be integrated by parts (it will change the coefficient of <f> 2 ). The remaining Lagrangian 
decouples into a dynamic part Coyn and an auxiliary (or algebraic part) Ca 



r v , r v 

'-Dyn t '-A 



(4.15) 



where 



C 



v 

Dyn 



Co C; 



2^2 2 
UJ 



L A 



4 

C 2 e" 2c 



-2a 



A 



v 



A v 

hYh v 



d r( j) v d r 4> v + e 



-2a ,V 2 



(4.16) 
(4.17) 



where A is defined in (4.12) 
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Equivalence with Regge- Wheeler 

We shall now demonstrate that in the 4d case the EOM derived from the Lagrangian 



( 4.16 ) turns out to be equivalent to the Regge- Wheeler equation || (see also |l2|) after 
a certain canonical transformation which essentially involves exchanging the field with its 
momentum. The Cunningham-Price-Moncrief master field ifcPM, m terms of which the 
EOM is usually written, is connected to our field <f) by 

C 

VCPM = — (toH? - C 2 C 2 d r (j> v ) , (4.18) 



2C 2 C 2 k v 



D=4 



where the auxiliary can be put to due to its algebraic EOM. Specializing ( |4.16 ) to 
D = 4 and Schwarzschild coordinates (|3.2|) the dynamic part of our action becomes 



S V Dyn = /dry^ + rVn- (4-19) 



with / = 1 — (the unimportant overall factor was changed to simplify the expressions). 

Before we perform the necessary canonical transformation we have to turn first to 
the Hamiltonian. The conjugate momentum is ir v = 5SDyn/8(d r (j) V ) = j^v $ V and the 
Hamiltonian is 

Now we perform a canonical transformation 

v 



^ = Z!_ ~v = _ r(j) v (4 21) 

r 

(in fact the new function <f) v is equal to (fcPM U P to a numerical factor). This transfor- 
mation is produced by the generating function F(cf) V ,4> v , r) = r <j) V 4> v . In general, for a 
generating function of the form 

F = r/(r)(t>4> (4.22) 
the canonical transformation is given by 

0=-, 7f = -# ; H = H + —<M>. (4.23) 
7] rj 

Performing the canonical transformation gives 7i = ^[A V (f> v2 — / _1 vr l/2 ) — 4> v tt v /r 
and then Legendre transforming back to a Lagrangian we find 

i = - - - (4.24) 

The corresponding EOM is 

^ V =(^- 3 -^W, (4.25) 



where □ = d r (fd r ) — f 1 u> 2 is the (r, t) Laplacian (or D'alembertian). This is precisely 
the Regge- Wheeler equation. 5 



5 Up to a change to Lorentzian signature which amounts to ui 2 
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These considerations can be easily generalized to the case of arbitrary dimension D and 
any background gauge. Starting with the Lagrangian (4.16), we turn to the Hamiltonian, 
perform a canonical transformation with the generating function 

F = e^ ^^ (4.26) 

according to ( |4.23| ) and transform back to the Lagrangian. The result is 

+ ^ e - 2a - Dc \Ae 2c uj 2 + 3(D - 2) 2 e 4a + (4C 2 - 2(D - 1)(D - 2))e 2a }^ 2 . (4.27) 
In the Schwarzschild coordinates this final form of the action encodes the master equation 

of § 5 

S v = / dT + { 4C 2 -2(Z>-l) (D - 2)/ _ 3^ + ^ ^ (4 28) 

where / is defined in (fT^). 
4.3 Scalar sector 

In this section we apply the same procedure to the scalar sector. In this sector there are 7 
fields: /if, /if,, h~ r , /if, /if,/i and h s . Their gauge variations are 

Shg = 2 (iu + a'e 2a £ Sr ) , Shg = £?' - 2a' £f + iu g , Sh s rr = 2 - b' g) , 

^ = 2c'e 2c e 5r > ^ = 2C S • (4.29) 

Altogether the number of fields and gauge functions is np = 7, «g = 3 and hence according 
to the general formula (^1|) we expect a single dynamic field (7 — 2*7 = 1) and 3 auxiliary 
fields. 

Indeed, in the terminology of Q the derivatively- gauged fields are hf r , hf r and /if, 
and the rest are algebraically gauged, so the general theorem says that the Lagrangian 
can be separated into an auxiliary part consisting of 3 fields H^ r , H^ r and iff (they are 
3 derivatively-gauged fields shifted by additional terms coming from a completion to a 
square) and a dynamical part consisting of a single gauge- invariant field cj) S . We define 



h s -2iojhf, (4.30) 



In terms of this field the Lagrangian can be written as 
S s = fe a ^ D - 2 ^ds£ s 



rS _ rS , /*S 
^ — l -Dyn ' J -A 



C Dyn = ^(d9^ s 5> 5 + /3v^^ 5 ' + 7^ 2 ) , (4.31) 
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where La is an auxiliary Lagrangian and (the expressions are organized by first expanding 
in powers of e 2c u 2 , each coefficient expanded again in powers of C 2 and a final expansion 
in powers of e 2a ). 



A = (D - 2) 2 eV + {D- 2){D - 3) C 2 (e 2a + 1) - 2(D - 2) e 

C 2 



'2a 



e 2c u 2 



+ ^-(D - 3) 2 \C 2 {e 2a + l) 2 - 4(2? - 2) e 2a 



(4.32) 



C2C2 



a 



(£>-2)(£>-3) 



^ e -3a-cf 2 p _ 2 ) [2C 2 - (£>-2)((£>- l)e 2a - D + 3) 



2c 2 

e a; 



+ 2Ci(D - 3)(e 2a + 1) - C 2 (-D -2){D- 3) [(3D - 11) e 4a - 2(D - 6) e 2a - D + 3 
+ 12(D - 2) 2 {D - 3) 2 e 2a (e 2a - 1)}, 



^ e _ 4a -2cf. 



2(D - 2) \2C 2 - (D - 2) ((£> - 1) e 2a - D + 3) J ^(D - 5) e 2a - D + 3j e 2c w 2 
+ (D - 3)|2C|(( J D - 5) e 4a - D + 3) 

- C 2 {D - 2) \2{D - 3){D - 5) e 6a - (3D 2 - 30D + 67) e 4a - 8(D - 3) e 2a + (D - 3) 2 



+ 4{D - 2) 2 {D - 3) 2 e 2a (e 2a - 1) 2 }V 



The auxiliary Lagrangian is given by 



C S A = ^H^MH S 



(4.33) 



with 



(H r s r \ 



(4.34) 



(the auxiliary fields H rr , Ht r and H r are defined in the Appendix |B|) and M is a 3 X 3 
Hermitian matrix given by 



M 



/ (g-2)(g-3) e 1a-Ab i{D _ 2) e a-36+c 



i°^±e a - 3b ~ c (e 2a + l) 



i(D - 2) e a - 3b+c u 
e 2a + 1) 



(-D-3) e a-3b-cf 2a 



-C 2 e- 2b 
-icj e~ 2fe 



zu; e 



-26 



Cof 1 e~ 2fe - 2c (2(D - 3) e- - <- x~ 



a 2a c 2c, ,2\ 



(4.35) 
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Equivalence with Zerilli 

We shall now show that just like in the vector case the standard Zerilli equation is obtained 
after an analogous canonical transformation. In D = 4 and Schwarzschild coordinates the 
master field of Zerilli (as given in p^]) is 

tq — r 



-7T 



3r + r{C 2 - 2) 
An action that yields Zerilli's equation is 

f^'z +Vz(r,io) *1 



Sz 



dr 



where 



and 



V z {r,u) 



A(r) 



^-2) 2 (£f + f?) + ^(C 2 -2 + ^) 



+ T 



A(r) = C 2 - 2 + 



3r 



(4.36) 



(4.37) 



(4.38) 



(4.39) 



Just like the case of the vector sector this Lagrangian is obtained from ours by per- 
forming a canonical transformation with the generating function 



A(r) 



(4.40) 



(and integrating by parts the term with one derivative). 

Generalizing to arbitrary dimension and background gauge the generating function is 



F 



where 



A (r)=ft + (C - 1) ' C - 2) (l-^) 



Then our usual procedure leads to the following action 



e a ds 



d r <f> s d r cj) s + V(r) 



is 2 



where 



-2c 



V(r) 



Af» 



Q{r)+e- Ia uj 2 , 



(4.41) 
(4.42) 

(4.43) 
(4.44) 



with 



Q(r) 
a 

7 



(D - 1)(D - 2) ^ 2a _ ^^ 2D 2 _ nD + 18 ) g 2a + (£> _ _ 6 ) J 5 

[D-l) 2 {D-2f 
16 



(e 2a - 1) 2 {(£> - 2) 2 e 2a - 2(D 2 -7D + 14)}. 
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Using Schwarzschild coordinates, we obtain the following action 




(4.45) 



which reproduces the equation obtained in |8[. 5 

5. Summary of results 

In this section we collect our main results. 

The background Schwarzschild metric, allowing for an arbitrary choice of the radial 
coordinate can be written as 

ds 2 {0) = -e 2a ^ dt 2 + e 2b ^ dr 2 + e 2c « dQ 2 D _ 2 , (5.1) 

thereby defining the background fields a, b and c. The perturbation field h^ u is defined 
through 

9w = 9$ + V ■ ( 5 - 2 ) 

The usual Lorentzian Schwarzschild metric ( |5.1|) has a Euclidean version ( |3.lD obtained 
simply by reversing the sign of gu- Since the metric is static the transformation between the 
two is quite straightforward and is effected by the analytic continuation t—*it. Whereas 
so far we used the Euclidean signature (for the sake of a certain mathematical uniformity) 
in this section we transform to the more physical Lorentzian signature. The effect of the 
analytic continuation t — » i t on our formulae is given by the following transformation 

u) — ► — iu 

h t -»• —iht ; h tr -> -ih tr ; h u -> -h a (5.3) 



The tensor sector. The tensor sector includes a single field and no gauge functions, so 
there is no issue of gauge choice or decoupling. 
The action in the tensor sector is given by 

S T = f e a+ ^ D - 2 ^dsC T , (5.4) 

where 

ds:=e b dr (5.5) 

is the Id "volume" element and e a+<yD ~ 2 ^ c ds = J \J—go dPx is the D-dimensional volume 
element reduced over time and the sphere. The Lagrangian density is given by 

C T = l [c 2 e~ 2a ^f* - d r 4> T d r 4> T * - V T {r) <f> 7 > T *] (5.6) 

where 

• (f) T = 4>^i(r) is the traceless and divergenceless part of hi; 
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V T (r) = C 2 (l)e 2c is the potential, and 
C 2 (l) = l(l + D-3). 



The vector sector. In this sector we count 3 fields and one gauge function, and therefore 
we have a single dynamic field (the master field) and a single auxiliary field. 
The action is given by 



S V = C 2 J e^ D ~^ ds [C v Dyn + C v A 



The auxiliary action is 



where 



and 



t~ A 



C2 H v jjVr* 



A v 

A v = C 2 (l)e-- - ,,r<: 



(5.7) 

(5.8) 
(5.9) 



C 2 := (Z-l)(Z + L>-2) . (5.10) 
The auxiliary field is an r-vector that can be expressed in terms of the fields hu V as follows 



H. 



A v h v r -iive- a {2e- c y/^ r -e- a d r )hl --C 2 e a - lc {2e- c ^r r -e- a d r )h v . (5.11) 



For the dynamic part we have two expressions. The first is the one that naturally 
emerges from the auxiliary /dynamic method and it is given by 



r v - 

*-Dyn 

where the dynamic field is 



Co O 



2^2 2 



-2a 



e- c (ojh v + 2i hi) 



(5.12) 
(5.13) 



(after dropping an overall (— i) factor that comes from the analytic continuation of (4.14)). 

The second expression encodes the standard Regge- Wheeler master equation (as gen- 
eralized by Kodama-Ishibashi) . It is 



S 



v 



e a ds 



oj 2 e~ 2a d> v 6 v * 



d r 4> v d r 4> v * -V v (r)4> v j> v * 



(5.14) 



where the potential in the vector sector is 



V (,•)- [l(D-2fe 2a + C 2 - {D 2Y \ r" 2 - 



(5.15) 



and the new dynamic field, which is essential the one canonically conjugate to <fi[ is 



iv _ ex P\— 



D - 2 c-b-a) 



A v 



0, 



r<f>Y 



(5.16) 
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The scalar sector. In this sector we find 7 fields and 3 gauge functions, and therefore we 
still have a single dynamic field (the master field) and this time there are 3 auxiliary fields. 
The action is given by 



S v = C 2 j e a+( - D -^ds [C s Dyn + C 



(5.17) 



The auxiliary sector includes the r-tensor H^ r as well as the two r-vectors H^ t and . 



Their action is a straightforward analytic continuation of (4.35,4.35) according to the rules 
(U). 

The first expression for the dynamic action C^ yn in terms of the master field <f)f t is 
quite involved and can be gotten by an analytic continuation of ( 4.31] ) . 

The second expression produces the standard Zerilli master equation (as generalized 
by Kodama-Ishibashi). It is 



S 



s 



e a ds 



2 — 2a iS iS* 



d r b d" 



is* 



V s (r) 



is is* 



where the potential in the scalar sector is 

rS 



V b (r) 



A 2 



Cl + aCl + (3C 2 +j 



together with 



(5.18) 



(5.19) 



(5.20) 



and the coefficients a, (3, 7 which depend only on D and e a ^ are given in ( 4.45 ). The new 
dynamic field (f> s is essentially the momentum canonically conjugate to (f>f t . 
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A. Derivation of the one-dimensional action 

The action for gravitational waves around a general background was given in ( |3.6| ). In this 
section we reduce it to obtain a ID action. First of all we need to separate the indices into 
sets according to their behavior w.r.t. rotations: t and r are denoted by a, /3, . . . while 
the angular coordinates are denoted by Latin letters i, j, . . . . In addition, we have to take 
into account the fact that covariant derivatives on the {D — 2)-dimensional sphere and in 
the ambient space differ in Christoffel symbols with at least one index in (r,t). In fact the 
only non- vanishing Christoffel symbols are of the form T^, r^ a , Tfj and We denote 
the covariant derivatives on a sphere by Vj. With indices separated in this way and after 
introducing covariant derivatives on a sphere the action gets the form 

S = J d D x^(A 3 + A 2 + A 1 + A + *_ (B 3 + B 2 + B 1 + B )) , (A.l) 

where the j4-terms and 5-terms come from the first and second terms in ( |3.6D , and sub- 
scripts highlight the number of (t, r) indices in the "parent" term. We have (where in what 
follows a comma means a usual partial derivative w.r.t. r, rather than a covariant one) 

A2 = V ih a pV 'ih a p + 2hi/3,afkp,a — 2Hip !a Kia,P — ^^ihaphip^a 

A = Vih jk Vih jk -2Vih jk Vjh ik + 4Tg% a V k hij - 4rgr^.% a ^, 
B2 = 2hpf3^ a hu t a + Vih aa Vihpp, 

Bo = V \hjjV 'ih kk . 

In these expressions the summation over the repeated indices is assumed to be done with 
the corresponding components of the inverse metric, which is diagonal; they are not written 
explicitly in order not to clutter the notation. So, for example, h aa should be understood 
as g^ a ^h a p, ha as g^^hij etc. Note that this action is valid for quadratic perturbations 
around any background of the form S m x X, where X and m are arbitrary. 

The explicit expressions for the inverse metric and the Christoffel symbols in the back- 
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ground ( jO| ) are 



g (0)tt = e -2 Q) g (0)rr = & -2h > ^(0)« = e -2c~^ (A 2) 

I*. = a', T[ t = -a'e 2 ^, T r rr = b> , I*. = </$?", rr. = -c'e 2 ( c ~%, (A.3) 

where ^ is a metric on a unit -D — 2-dimensional sphere. 

Next we separate the a, j3 indices according to a — > (t,r), express ft's in terms of /i's 
according to the definition ( |3.4|) , expand the fields as in (^) and calculate angular integrals 
with the formulae of Appendix (|C|). In this way one arrives at the action described in 
subsections [0|-[0[ 



B. Auxiliary fields in the scalar sector 

In this appendix we give the precise definitions of the auxiliary fields in the scalar sector. 
As mentioned in section |4.3| there are 3 such fields that come from a completion to a square 
of fields h^ r , hf r and hf . We called these fields H^ r , H^ r and . They are invariant w.r.t. 
the gauge transformations ( f4,29j ) and appear in the Lagrangian of the scalar sector in an 
algebraic way (without derivatives). 
The first of these fields is given by 

H^r := e ~ 3a ~ c+2b (aihf r + a 2 hf t + a 3 hf + a^h s + a 5 h s + a 6 /if/ + a 7 hf + a s h s ' + a 9 h s '^j 

(B.l) 

where the expressions that follow for the coefficients are organized (as in section ^.3|) by 
first expanding in powers of e 2c uj 2 , each coefficient expanded again in powers of C 2 and a 
final expansion in powers of e 2a . 

at = e 2 ^ A s 



a 2 = ir~ = - — {D - 2 ) 



2uj 2 
C2, 



(D - 5)e 2a -D + 3 



e 2c uj 2 



(D - 3)(e 2a -1){C 2 {D- 5)e 2a + D 



2d 



4(D-2)(D-3)e 



2a 



«4 = -e 



-2c 



2(D - 2f (D - 5)e 2a -D + 3 e 4c a; 4 — (D — 2)(D — 3)e 2c L0 2 x 



„2c .2, 



I C 2 {D - 13) - A{D - 2){D - 5) e 4a + 2 C 2 + 2{D - 2){D - 3) e 2a - C 2 (D - 3)} 



+ 4C 2 (D - 3) 2 C 2 e 6a 
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«5 = C 2 C 2 ^-^e 4a - 2c {2(D - 2)e 2c u? + C 2 (D - 3)e 2a }; 



a 6 = ^ = -^e a - b+c l 2(D - 2)e 2 V + (D - 3) 



C 2 (e 2a + 1) - 2(£> - 2)e 



2« 



a 8 = ~^e a - b - c l 2{D - 2) W + (D - 2){D - 3) 



C 2 (3e 2a + 1) - 4(L> - 2)e 



2« 



2c 2 

e cu 



+ C 2 C 2 (D-3) 2 (e 4a + e 2 ^ 



«9 



C 2 C 2 D - 3 



2 D 



e 3a-6-c| 2 (D - 2)e 2c u; 2 + C 2 (L> - 3)(e 2a + 1)}. 



The second auxiliary field is given by 

a - c (ft /l& + ft ^ + ft /if + /?4 /i 5 + ft ^ + ft /4' + 07 hf + ft h S ' + ft h S ') , 

(B.8) 

where the coefficients are 

ft = A s 



ft = ^e- a+b+c u;l 2(D - 2)(e 2a - 1) [(D - 2){D - 3) - C 2 



e 2c o; 2 



- (D — 3){c 2 2 (e 2a + 1) - C 2 {D - 2) [2(D - 4)e 4a — (D — 7)e 2a -D + 3 
+ 4(L>-2) 2 (D-3)(e 4a - e 2a 



4 



ft = -^e- a+b - c |4(D - 2)e 4 V 

+ 2(D - 3) [c 2 (e 2a - 1) + (£> - 2)(2e 4a + (£> - 7)e 2a -D + 3) 
+ {D- 3) 3 (e 2a - 1) rC 2 (e 2a + l) 2 - 4{D - 2)e 2a 



e 2c u 2 
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p 4 = - l -e- a+b - c uj^2(D - 2) 2 e 4c w 4 



e 2c u 2 



+ {D — 2){D - 3){c 2 (3e 2a + 1) + (D - 2) 2e 4a + (D - 9)e 2a - D + 3] } 

+ (jP ~ 3)2 {2C 2 2 (e 4a -e 2a ) 
+ C 2 \(D + 5)(£> - 2)e 6a + (D — 2)(D — 13)e 4a - (D - 1)(D - 2)e 2a — (D — 2)(D - 3) 

- 4{D - 2) 2 e 2a (2e 4a + (£> - 5)e 2a - D + 3)} | ; 



ft = -yC 2 GD - 3)e^" c 2e 2c u, 2 + §-| 



C 2 (e 2a + l) + 4(D-2)(e 4a -e 



4a 2a > 



/5 6 = 2) J 2(D - 2)e 4c a; 4 + ID - 3) C 2 (3e 2a + 1) - 4(D - 2)e 

4a; 1 



2a 



e 2c u 2 



(3 7 = ^-(D - 3)1 2(D - 2)(e 2a - l)e 2c uj 2 - (D - 3) 



C 2 ( e 2a - l) 2 -4(£>-2)e 



2« 



/% = - 2)(D - 3)(e 2a - 1) j (£> - 2)e 2c cu 2 + 



C 2 (3e 2a + 1) -4(D-2)e 



2a 



/ gb = ^C 2 (£)-3) 2 e 4a a;. 



The third auxiliary field is 

#S = C 2 e- fl - c ( 7l ^ + 72 /4 + 73/tf + 74 /i s + J 5 ~h s + 76/tf/ + l7 hf + l8 h s ' + l9 h s '' 

(B.17) 

where the coefficients are 

71 = A S 
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273 J- - 

72 = 2^7 = ~r 



-a+b+c 



2(D -2) C 2 - 2(D - 2)e 



2a 



e 2c co 2 



+ C 2 (D- 3){c 2 (e 2a + !)-(£)- 2)e 2a [(£> - 3)e 2fl - L> + 5] } j ; 



74 = ~-e 



-a+b—c 



2{D - 2) 2 e Ac uj A 



+ (D — 2){D - 3){c 2 (3e 2a + 1) - (D - 2)e 2a (D + l)e 2a - D + 3 }e : 
+ C 2 {D - 3) 2 [C 2 (e 2a + !) + (£>- 2)(e 4a - 3e 2a ) 



2c a; 2 



75 = ~{D - 3) e a+6 - c \ 2 \C 2 - 4(D - 2)e 2a 



e 2c u 2 



+ C 2 



D-3 
D-2 



C 2 (e 2a + 1) - 4(D - 2)e 



2a 



76 



Z7 7 



D-2 



2c 



2tu (D - 3) (e 2a - 1) 



[2{D - 2)e 2c uj 2 - C 2 (D - 3)(e 2a - 1)}; 



79 



C2 
2 



(D - 3) e 2a {2(D - 2)e 2c uj 2 + C 2 (D - 3)}; 



C. Multidimensional spherical harmonics 
General theory 

We consider tensor spherical harmonics of various ranks on a d— 1 dimensional sphere, where 
in our application d = D — 1. We give here the properties that we need, while furthermore 



an explicit construction of the harmonics can be found, for example, in [16]. Spherical 
harmonics transform in an irreducible representation of SO(d). These representations are 
labeled by Dynkin indices [Zi.../ n ] where n = [d/2] is the rank of the group SO(d) and 
each function within the representation is further labeled by a "weight vector" [mi...m n ]. 
Henceforth we shall denote the multi-index m\...m n collectively as m. 
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In this context it is necessary to distinguish two types of tensor types. The first 
tensor-type is with respect to local diffeomorpshims under GL(d) and it will be denoted 
as a fi-scalar-field or ri-vector-field etc. The second tensor type is with respect to the 
SO(d) global isometries: we will refer to [Z, 0, ...0] as an SO(d) -scalar, to [Z, 1, ...0] as an 
50(c£)-vector and to [1,2, ...0] as an SO(<Z)-tensor. 

First, there are (Q) scalar spherical harmonics Y lm that transform under rotations of a 
sphere in the SO{d)-sc8lax representation [1,0, ...0]. These functions are all eigenfunctions 
of the Laplace-Beltrami operator □ = VjV* 

UY s lm = -C 2 Y s lm , (C.l) 

where C 2 ■= 1(1 + d — 2) is the eigenvalue of the quadratic Casimir operator in the SO(d)- 
scalar representation [7,0, ...0]. It is a generalization of the well known formula for d = 3 
spherical harmonics OY = —1(1 + 1)Y. 

Then there are (£1) vector spherical harmonics, which belong to two families according 
to the Hodge decomposition. The first family consists of O-vector-fields of the following 
gradient form 

Yi S im = ^iY lm . (C2) 

These are SO (d)-scalars which explains the superscript S (they belong to the same type 
of SO(d) representation [Z, 0, ...0] as the fi-scalar-fields) . In view of (CA) they obey 

V%Y lr £ = -C-i Y lm , (C.3) 

The functions in the second family of (O) vector harmonics YY, belong to the SO(d)- 
vector representation [Z, 1, ...0] and are diver genseless: 

V^L J = • (C.4) 

The action of the Laplace-Beltrami operator on Yf can be calculated from the definition: 

a Y* lm = -(C 2 -d+2)YS lm . (C.5) 

whereas for Y^ the general theory presented in gives 

UY^ lm = -(C 2 -l)Y^ lm . (C.6) 

Then there are (Q) symmetric tensor spherical harmonics, which come in 4 families. 
First of all, as gij and Vj are invariant tensors, it is obvious that the f2-tensor valued 
functions 9ijY s lm and VjVjYf m transform in the SO(d) -scalar representation [Z, 0, ...0]. 
However they are not orthogonal w.r.t. to the inner product 

< Yi\Y 2 >= J Y 1Aj Y* j dn . (C.7) 
S d-i 

In order to orthogonalize this system of functions we define two (O) tensor spherical har- 
monics by 



Y 



C 2 

' ~ 9ij Y i m , lm := (VjVj + - — ^}Y lm . (C.8) 



ij, lm • 
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Then there are spherical harmonics of the gradient form 



Y: 



v 



1 



l i, lm) ' 



(C.9) 



which belong to the vector representation [Z, 1, 0...0] and are traceless because of (|C.3|) . 
Finally there is the family of (Q) tensor spherical harmonics lm which belong to the 
symmetric traceless tensor representation [Z, 2, 0...0] and is divergenceless. The action of 
the Laplace-Beltrami operator on Y s , Y s and Y v can be calculated directly from their 
definitions: 



□X- 

^~^~ij, lm 



ij, lm < ~~'2 1 ij t l m ? 



-{C 2 -2d + 2)Y i s j 



' ij, lm (C2 

and the action on Y T is given by (see [|16f) 



d - DY 



ij, lm s 
V 

ij, lm ' 



□Y- T , 

ij, lm 



(C 2 -2)Y i 



ij, lm " 
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Divergences of tensor harmonics are given by 
VjY Vij = _C 2 -(d-l) 



d-2 



(C 2 -(d-l))Y s i 



d-l 

Y Vj , V,Y Tij = . 



(C.12) 



Orthogonality and normalization 

The spherical harmonics constructed above are orthogonal to each other. For harmonics 
belonging to different irreps it follows from the usual theorems about unitary represen- 
tations, and Y$ and Y$ are orthogonal by construction. Now we calculate the norms of 
various spherical harmonics, where in fact the first line in each group is just the definition 
of the appropriate normalization. 
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Spherical harmonics on S 2 

The case of spherical harmonics on S 2 is special. The rank of SO (3) is 1, so "weight vectors" 
consist of a single number, and therefore all harmonics, no matter if they are scalar, vector 
etc, belong to the same representations of SO (3). Scalar spherical harmonics are labelled 
by single Dynkin index I and a single weight m: Yi m . Vector spherical harmonics can be 
built out of scalars with the help of invariant tensors. There are 2 independent invariant 
tensors: the covariant derivative on the sphere Vj and the Levi-Civita tensor e^. One can 
build 2 families of vector harmonics: 

= V 4 Y W , Y}~2 = e/V,lW . (C.16) 

The sign in the superscript denotes the spatial parity of the harmonic with I = 0. Higher- 
dimensional spherical harmonics Yf and YY are generalizations of and ' corre- 
spondingly. 

One can continue the process of covariant differentiation and multiplication by to 
build tensor spherical harmonics of higher ranks. For example, one gets four families of 
rank 2 tensor spherical harmonics: 

YtSij = *?ViVMm , ^% = eiS n VkV n Y lm , (C.17) 

8 families at rank 3 etc. For our purposes we need vector and symmetric rank 2 tensor 
harmonics. Out of four rank 2 tensor harmonics one can build three symmetric ones: 

I&S = > = e ( h£VkVnY lm , Y^\ = e ( f V^V^ , (C.18) 

where the first two have even parity while the last one is odd. 

The spherical harmonics Y^^ and Y^^j in ( C.18 ) are not orthogonal w.r.t. the inner 
product ( C.7 ) (harmonics of the opposite parity are obviously orthogonal). To orthogonal- 
ize them we notice first that 

+ y i +2) = ~9isKl + VY ■ (C19) 
Therefore the orthogonal linear combinations are similar to 

Y {+) -a- Y, 

Y lmlj = ViVjlj m - \diPY lm = ViV 3 Y lrn + l S^±}lg ij Y lm . (C.20) 

In higher dimensions the harmonics Y^\ Y^ and y}- "* become Kf, Kf and YY . There 
is no analog of Y? on S 2 . 

The actions of Laplace-Beltrami operators on various spherical harmonics are directly 
obtained from the higher-dimensional formulas ( C.10| ), their divergences from ( C.12 ) and 
the normalizations from ( p. 13 - C.15[) . 
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